Mounding Instability and Incoherernt Surface Kinetics 
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Mounding instability in a conserved growth from vapor is analysed within the framework of adatom 
kinetics on the growing surface. The analysis shows that depending on the local structure on 
the surface, kinetics of adatoms may vary, leading to disjoint regions in the sense of a continuum 
description. This is manifested particularly under the conditions of instability. Mounds grow on 
these disjoint regions and their lateral growth is governed by the flux of adatoms hopping across the 
steps in the downward direction. Asymptotically lnt dependence is expected in 1+1- dimensions. 
Simulation results confirm the prediction. Growth in 2+1- dimensions is also discussed. 
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Mounding instability was experimentally observed and 
proposed by Johnson et. al fm during growth of GaAs 
on (001) GaAs substrate. Initially, activation difference 
(Schwoebcl - Ehrlich (SE) barrier ) between adatoms 
hopping on the plane and the one crossing the step edge 
was considered responsible ||. Later it was shown that 
edge diffusion can also lead to similar effect j|. One 
of the issues related to growth of mounds has been the 
temporal dependence of mound growth. Based on vari- 
ous forms of continuum equations the lateral growth is 
expected to have a time dependence ~ t s where s takes 
values from 0.0 to 1/4 ||. Similarly, the width of the 
interface is predicted to follow the power law t" with j3 
varying from 1/3 onwards S . All these predictions are 
based on the assupmtion that the underlying conserved 
growth equation describing non equilibrium growth, is 
valid over the entire substrate. In the following we show 
that under the conditions of instability and low temper- 
atures, this assumption fails. Consequently the growth 
of mounds is governed laterally by the adatom kinetics 
across the mound boundaries. We show this by devel- 
oping growth equation over a stepped surface in 1+1 di- 
mensions using kinetics of adatoms and steps. This helps 
establish the correspondence between kinetic processes 
and terms in the growth equation. We assume that only 
mechanism of relaxation is by diffusion of adatoms. This 
allows identification of process and corresponding term 
uniquely. These assumptions are expected to be valid at 
low temperature, where evaporation is negligible. Once 
the kinetic processes leading to various terms in growth 
equation are identified, presence or absence of such terms 
in various regions on the surface can be predicted. This 
allows us to classify different regions on the surface ac- 
cording to the growth equation followed there. 

Consider growth on a one dimensional substrate. Fig.l 
shows the stepped region under consideration. Growth 
proceeds through randomly falling adatoms on the sur- 
face that relax by diffusing on the stepped terraces. 
Adatoms with zero nearest neighbors (nn) are mobile 
while those with more than zero nn will have negligi- 
ble mobility. Further, desorption and dissociation from 



the steps is also negligible at low temperature. Under 
the conserved growth conditions it is possible to write 
formally the growth equation in the form dth(x., t) = 
V • j(x, t) + F, where, F is incident flux, h(x, t) is height 
function and j(x, t) is particle current. An uphill current 
on a tilted substrate indicates instability while downhill 
indicates stable Edward- Wilkinson (EW) || type growth 
. Let Id be the average length travelled by an adatom 
before getting attached to another adatom or step. The 
density of steps can be expressed as 1 .j^L . Let Pa and 
Pb be the relative probabilities for hopping across the 
sites A and B in Fig. 1. By considering current due to 
the downward hops and that due to the in-plane hops 
seperately, one can show that the resultant nonequilib- 
rium current is given by flq] 
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Where n denotes +ve x direction. Presence of in the 
denominator accounts for the nucleation effect on larger 
terraces. In this expression, local terrace width is (l^ 1 + 
a~ |m|) _1 . However due to the relative velocity between 
two adjacent terraces, the local terrace width changes. 
The velocity difference will be proportional to the ^gj^ • 
Including this dynamical effect, the expression for the 
current becomes, 
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Next, we argue that every downward hop introduces 
height-height correlation, hence will give rise to all the 
stabilizing terms in a growth equation. Under the tilt 
independent current conditions, the lowest of such terms 
is K§£%£. Thus the current on the stepped surface will 
be, 
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For small slopes, above current generates growth equa- 
tion in the moving frame with average growth rate, 

dh(x,t) = F(P B - P A )l d d 2 h 
dt 2 dx 2 

Fli d 2 fdh\ 2 B 4 h , , 

+ -T^UJ - fc ^ + ^' t) (4) 

Where, ?y(x, i) is the Gaussian noise in the deposition 
with the property, < rj(x' ', t')rj(x, t) >= 8{x' — x)S(t' — t). 
For Pa = Pb, current is tilt free, corresponding equation 
has the Lai- Das Sarma - Villain HH form. 

Now consider a top terrace. By definition, only tilt 
independent terms will contribute. Also since steps are 
absent the only term that contributes is, fcfjr due to the 
downward hops at the edges. For a base terrace, only in- 
plane hops are possible, hence none of the above terms 
contribute. This shows that these regions offer restricted 
kinetics, thereby changing the growth equation locally. 
If these regions are smaller than I4, they will act only 
as the regions of discontinuity. Villain JhJ encountered 
this difference in connection with simulations using Zeno 
equations. However so far it has been assumed that a 
given growth equation is valid over the complete sub- 
strate and such a distinction is irrelevent under coarsen- 
ing. Based on above argument we propose that stepped, 
top and base terraces are distinct regions where different 
growth equation applies. Thus, in growth from vapor at 
low temperature, scalability breaks down. However, in 
simulations and experiments, kinetically rough surfaces 
are observed. This is possible provided steps develop 
in base and top regions at a rate comparable with the 
growth of correlation length defined by the equation over 
steps, shifting these regions. It may be noted that since 
downward hops allow addition of material to lower layers, 
average shifting of top or base region is possible only if 
downward hops are present. An infinite SE barrier leads 
to such an immobilization of top and bases Jll| . This 
is the limiting example displaying the effect of three in- 
dependent regions on growth. If average shift of top or 
base regions lags behind the correlation growth, mound- 
ing instability appears. The lateral growth of the mounds 
being decided by the lateral shift of base regions. It has 
been established that stability of the growth can be quan- 
tified in terms of tilt dependent current j t . For j t > 
(uphill), unstable growth while for j t < 0, stable growth 
is obtained, where jt is measureable in a simulation by 
properly adjusting the boundary conditions 0. Thus, 
for uphill current, base or top regions lag behind the cor- 
relations while for downhill it shifts at least as fast as 
correlations. 

Consider unstable growth in 1+1 dimensions. Fig. 2 
shows the well developed mounds. Note the deep ridges 



formed due to high step heights of the steps forming the 
ridges. The model used for this growth will be described 
later. It suffices to know that its a growth with finite 
diffusion of adatoms and finite SE barrier. We estimate 
the growth rate by appealing to the diffusional kinetics 
of atoms. The growth proceeds by expansion of a larger 
mound at the cost of smaller one pi . Thus the ridge pro- 
ceeds in one direction. A smaller mound generally poses 
a smaller angle with respect to the substrate. Thus, rel- 
atively longer terraces are present on this mound. The 
diffusional addition to the ridges is mainly from these 
terraces, resulting in to the shift of ridge in the direction 
of smaller mound. Thus, we assume that adatoms are 
added from the smaller mound, diffusionally. The diffu- 
sional rate of displacement is dl = D 1 s ^ 2 t^ 1 / 2 dt on a plane 
surface in time dt. However, for a ridge to move laterally, 
it must be filled at least up to first step height. For sharp 
ridges as in Fig. 2, the step height of ridge may be taken to 
be ~ w, the rms height fluctuation (width). Hence, the 
displacement for a ridge will be dl r — pD\^' 2 t~ 1 l 2 aw~ 1 . 
Where, p is relative fraction of adatoms crossing the step 
edge and a is lattice constant. For w ~ i 1 / 2 , the growth 
of mounds is proportional to Int. 

We verify the Int dependence for a I+I dimensional 
model that mimics the growth at low temperature. In 
this model, on a one dimensional substrate, adatoms are 
rained randomly. An atom with one or more nn is incor- 
porated in the crystal. An adatom with zero nn is allowed 
to hop n number of times at the most. If it acquires a 
nn, then no further hops are allowed. If number of hops 
are exhausted, it is incorporated at the final site after n 
hops. A parameter p is introduced, such that for p > 0.5 
hopping across a step in the downward direction is dif- 
ficult, p = 1 is the case of infinite SE barrier. We have 
measured < hihj > correlations for various values of p 
and used the first zero crossing as the measure of the size 
of the mound. In Fig. 3, plot of mound size Vs. time on 
a semi log scale clearly shows that for p > 0.5 i.e. for pos- 
itive SE barrier, the mounds growth is Int. Also shown 
is the case for p = 0.5. We plot length corresponding to 
the first maximum in height-height correlations for this 
case. The curve on semi-log plot is exponential showing 
a power law dependence. Correlation length £ ~ t 1 / 4 in 
this case. In fact it can be shown || that correspond- 
ing equation describes Das Sarma - Tamborenea (DT) 
|12| model to which the tilt independent growth equa- 
tion reduces for large slopes. We find that for p < 0.5, 
asymptotically, EW growth is recovered. Thus, the base 
and top regions move at least in phase with the £ to 
provide rough surface. In the present model dissociation 
from steps is not included so that the detailed balance is 
not followed. If this is included, and the current is still 
uphill then Int dependence continues for growth in 1+1- 
dimension. 

Above arguments are true in any dimension. In 2+1- 
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dimensions, mound formation is observed experimentally 
as well as in simulations p]Jl3]|. Various predictions are 
referred in the introduction above regarding the time evo- 
lution of the mounds. The Int dependence in 1+1- di- 
mensions is the upper limit for lateral development of 
the mounds in 2+1- dimensions. This is so because, 
a given mound is surrounded by four or more mounds. 
Probability that such a mound happens to be the small- 
est amongst the surrounding ones including itself is very 
small. A given mound may be reduced in one direction, 
but it may increase in other direction owing to a smaller 
mound there. Thus, instead of consumption, shift of 
mounds is more likely on a two dimensional substrate. 
In order to find the time dependence of mound growth 
in 2+1- dimensions, we have used same model described 
above, except that the rules apply in two directions on 
a square lattice. In addition, we have included edge dif- 
fusion with no edge barriers. It is observed that edge 
diffusion suffices to induce uphill current so that even if 
the diffusion of single adatoms is unbiased, mound for- 
mation is observed. In the absence of edge diffusion but 
with unbiased single adatom diffusion, EW type growth 
is obtained |Q. Noise reduction technique [|l|is em- 
ployed with reduction factor of 5. The growth of mound 
size is monitored in the same way as for the 1+1- dimen- 
sions, using zero crossing for the correlations < hihj >. 
Fig. 4 shows the plot of mound size as a function of time 
on semi-log plot. Clearly, after an initial growth like 
Int, the curve tends to saturation, confirming the slower 
growth rate. By varying parameter p, a condition close 
to tilt independent current is obtained. The growth in 
that case follows, t 1 ^ power law. From the arguments 
leading to Eq.4 , in 2+1- dimensions, we find that asym- 
metric term will be ineffective if step edge tension is lower 
so that steps morphology is wavy or fingered. This is so 
because, the terrace size can be reduced by step move- 
ments in the orthogonal directions as well. Thus only 
V 4 ft term contributes, leading to (3 — 1/4 and z = 4 in 
2+1- dimensions. Clearly, this observation suggests that 
in experimental growth, if SE barrier is very small(but 
nonzero), at low temperature growth rate of mounds can 
be i 1 / 4 in the transient region. If the edge tension is high 
so that steps are straight and less wavy, asymmetric term 
can contribute with (3 and z, characteristics of a Lai- Das 
Sarma like equation Q in the transition region. 

In conclusion, we have shown that growth from vapor 
on surface proceeds via in principle a heterogeneous dy- 
namics. The stepped, base and top regions on the surface 
allow different growth dynamics. As a result the spatial 
scalability breaks down. The effect is distinctly observ- 
able for unstable growth leading to mound formation. 
The kinetics across the mounds suggest a Int dependence 
in 1+1- dimensions which is verifiable in a suitable model. 
A slower growth is predicted in 2+1- dimensions which 
is also observed in a model simulation. 
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FIG. 1. A typical step structure formed during growth 
along positive slope, v and v' are velocities of the steps . 
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FIG. 2. Morphology of the surface in 1+1- dimensions for 
an unstable growth after 10 6 layers. Parameter p is 0.6. 
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FIG. 4. Shows time evolution of lateral growth in 2+1- 
dimensions. The values of parameter p are 0.35, 0.6, and 0.7 
respectively for the curves from top to bottom in the figure. 
The substrate size is 300 X 300 for the simulation. 
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FIG. 3. Shows time evolution of lateral growth in 1+1- 
dimensions. The values of parameter p are 0.5, 0.6, 0.7 and 0.8 
respectively for the curves from top to bottom in the figure. 
The substate size is L = 10000. 
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